Superradiance of low density Frenkel excitons 
in a crystal slab of three-level atoms: Quantum interference effect 
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We systematically study the fluorescence of low density Frenkel excitons in a crystal slab con- 
taining Nt V-type three-level atoms. Based on symmetric quasi-spin realization of SU(3) in large N 
limit, the two-mode exciton operators are invoked to depict various collective excitations of the col- 
lection of these V-type atoms starting from their ground state. By making use of the rotating wave 
approximation, the light intensity of radiation for the single lattice layer is investigated in detail. 
As a quantum coherence effect, the quantum beat phenomenon is discussed in detail for different 
initial excitonic states. We also test the above results analytically without the consideration of the 
rotating wave approximation and the self-interaction of radiance field is also included. 
PACS numbers: 42.50 Fx, 71.35-y 



I. INTRODUCTION 

It is well known that the fluorescence of an exciton ex- 
hibits super-radiant character due to the appearance of 
macroscopic transition dipole moment of the exciton . 
However this collective feature of exciton radiance de- 
pends upon dimensionality of crystal. In infinite bulk 
crystals, the exciton does not radiate because this exci- 
ton is dressed with a photon which has the same wave 
vector due to the transitional symmetry of the total sys- 
tem. As a result, a stable polariton is formed 1 . In the 
case of lower dimensional systems, one can show that the 
exciton decays superradiantly due to breakdown of the 
translational symmetry. For example the decay rates of 
exciton are of the order (A/a)7 for one-dimensional (ID) 
crystals and (A/a) 2 7 for 2D crystals 2 , with 7 being the 
radiative decay rate for an isolated atom, a the lattice 
constant, and A the light wavelength. 

The enhanced factor (\/a) D in the radiative width of 
exciton is now regarded commonly as one of evidences of 
"superradiance" 10 , which has been demonstrated in ex- 
periments for Frenkel excitons 13 and Wannier excitons 14 . 
Since the enhanced radiant effect can also appear in 
the usual superradiance of an ensemble of atoms 15 ' 16 , 
how can we distinguish the different features between 
the fluorescence of the excitons and that of the atomic 
ensemble? It is known that, in the fully population- 
inverted atom systems, the atoms radiate independently 
with each other in the initial stage. The back-action of 
the emitted photon to the atoms results in the corre- 
lation among atoms. Consequently, the atoms become 
cooperative and thus the fluorescence from the atomic 
ensemble will show different statistical properties in the 
initial and final stages. On the other hand, the exciton 
fluorescence exhibits identical statistical character during 
the whole process 17 . Physically, this is because the ini- 
tial dipole moments of the atoms are spatially random in 
an atomic ensemble, but in a semiconductor crystal, the 
dipole moment of exciton presents a macroscopic effect 
even at initial moment. So the optical properties of multi- 



dimensional quantum-confined semiconductor structures 
(MQCS), such as quantum wells, quantum wires and 
quantum dots, have their own specific features in physical 
processes. 

The optical properties of the MQCS in a semiconduc- 
tor microcavity (SMC) have attracted more and more 
attention in the past decades 18 . The SMC with high- 
reflectivity dielectric mirrors leads to the realization of 
the strong coupling between radiation and matter. More- 
over, since the optical mode structure of the SMC may al- 
ter around the MQCS, many new phenomena, such as tai- 
loring the spontaneous-radiation rate and pattern 19-21 , 
the coupled exciton-photon mode splitting in a SMC 22 ' 23 , 
have been demonstrated. The resonant interaction be- 
tween excitons and a single-mode cavity field and the cor- 
responding detuning effect were further investigated 24 ' 25 . 

Most of these former works mentioned above dealed 
only with the two-level lattice atom case, however the 
three-level atom case may be very useful to implement 
quantum information encoding and processing 26-30 . 
Over the last few years, the cavity QED with the col- 
lective excitations of ensembles of three-level atoms has 
attracted much attention for quantum computing im- 
plementations. In this case, many atoms are entangled 
through their interaction with the common cavity field. 
To maintain quantum coherence in this quantum infor- 
mation processing 31 , it is important to reach the so-called 
strong coupling regime where the single-photon coherent 
coupling go 3> 7> 7cav, the atomic and cavity dissipa- 
tion (decoherence) rates respectively. It is the symmet- 
ric collective excitation that can reach the strong cou- 
pling regime without requiring a high finesse cavity as 
go cx \/]V, with the total number of atoms N and here we 
need to consider the quantum decoherence induced by the 
spatial inhomogeneousness of couplings 32,33 . In fact, the 
phenomena of superfluorescence or superradiance 34-36 
constitute another example of collective state dynam- 
ics. Recent experimental success clearly demonstrates 
the power of such an atomic ensemble based on the sys- 
tem for entangling macroscopic objects 37 . 
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In the present paper, we study the fluorescence of low 
density excitons in a crystal slab containing V-type three- 
level atoms. The purpose of this paper is to investigate 
the quantum interference effect 38 in the time evolution 
of light intensity. In sec. II, considering that the exis- 
tence of the two-level atomic exciton is mathematically 
associated with the infinite dimensional reducible repre- 
sentation of SU(2) Lie algebra 39 , we can rationally de- 
fine the exciton operators for the three-level case associ- 
ated with SU(3) algebra, which contains various SU(2) 
subalgebras. With this conception, both the free part 
and the interaction part of Hamiltonian can be written 
down in terms of the introduced two-mode exciton op- 
erators, which can be described by bosonic operators in 
the low-density limit. In sec. Ill, with rotating wave 
approximation (RWA) for the interaction (e/mc)P ■ A, 
the coupled equations for the exciton-light field system 
are solved with the help of Weisskopf-Wigner approach 
(WWA) 40,41 . In sec. IV, the explicit expressions of the 
electric field operators are derived for the monolayer case. 
The light intensity as well as the first- and second-order 
degree of quantum coherence are calculated to show cer- 
tain new features of exciton fluorescence in a three-level 
crystal slab. We discuss the phenomenon of quantum 
beat in the time evolution of the light intensity for vari- 
ous initial exciton states. In sec. V, we consider the roles 
of both the non-RWA terms and the self- interaction term 
of photon. This consideration avoids "unphysical" roots 
of the characteristic equation when the non-perturbation 
approach given in Ref. 17 is used. The light intensity is 
calculated to compare with the results in sec. IV. The 
conclusions are presented in Sec. VI with some remarks. 



II. TWO MODE EXCITON SYSTEM WITH SU(3) 
STRUCTURE 



We consider a plane crystal slab with a simple cubic 
structure which contains a stack of N identical layers. 
V-type three-level atoms, as shown in Fig. 1, occupy Nt 
lattice sites, where Nt = N^N and Nl is the total lattice 
sites within each layer. The wave vectors of the excitons 
and light fields are all assumed perpendicular to the slab. 
We restrict ourselves to investigate only the low density 
excitons region. 




(b) 




FIG. 1. (a) Energy structure of a V-type three-level atom, 
fliis the transition frequency of \g) <-> |ei), and f^2 is that 
of \g) <-> |e2). The direct transition between the two upper 
states is forbidden, (b) A plane crystal slab with simple cu- 
bic structure containing a stack of N identical layers with Nl 
sites within each layer. 

The interaction — p- A between the radiation field and 

mc 

the multi-atom system is written in the second quantiza- 
tion form 



(fj\ 



Jqla 



E^) |e 2 ) y (g\ 



iqla 



H.c. 



(1) 



where a q and a q 



q are the annihilation and creation op- 
erators of the photon with wave vector q along z, re- 
spectively, j and / denote the jth lattice site in the Zth 

layer. g x {q) = yf^fctk, g 2 (q) = y/^fc<h are the ef- 
fective atom-photon coupling constants of \g) «-> \e\) and 
\g) *-> \e2) atomic transitions with transition frequencies 
fix and CI2, respectively. The corresponding transition 
dipole moments are represented by di and d2. We as- 
sume that di and d 2 are parallel and lying in the slab 
plane. 

Introduce the collective operators for the Ith layer as 17 



l(0t = 



m = 




It is easy to prove that, two sets of the operators 
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and 
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(3) 



(4) 



just generate two SU(2) algebras not commuting with 
each other. This means the four collective operators 



Ei, E% , Fi' and F^' do not span a product alge- 
bra SU(2)<£>SU(2). A straightforward calculation in Ap- 
pendix A checks that they satisfy SU(3) algebra. Actu- 
ally, the above four collective operators define a spinor 
realization of SU(3) of Nl + 1 dimensions. Furthermore, 
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the unique number Nl, atom number, determines the 
dimensions Nl + 1 of representations of the two SU(2) 
subalgebras. Since we understood the single-mode ex- 
citon in terms of the large Nl limit of representations 
of SU(2), it is easy to prove for SU(3) case that, in a 
very large Nl and the low density excitation region, the 
above defined collective operators and B^> become 
the bosonic ones obeying the commutation relation 



^(0 A (l ' n 



Su 



(5) 



where the ideal bosonic approximation is equivalent to 
the neglect of the phase space filling effect and the 
exciton-exciton interaction 42 . 

In terms of these collective operators, the two-mode 
Frenkel exciton operators with wave vectors k = 
(to = -±(JV - 1), -±(7V - 3), l(N- 1)), are just the 
discrete Fourier transformations for them 

Vn^ 



(6) 



In fact their conjugates A\ and Bl. are just the generators 
for the quasi-spin wave states 



\ A k) = A\\g 1 ,g 2 ,...,g NT ), 
\ B k) = B k \g u g 2 ,...,g NT ). 



(7) 



Since the operators Ak and B k commute with each other 
for very large Nl and low excitation, they form a two 
independent boson system. Then we obtain the interac- 
tion Hamiltonian for the the two-mode Frenkel exciton 
system coupled to a quantized electromagnetic field 



H x =h^Gi(q)0(k + q) 

q,k 

+ hJ2G2(q)0(k + q) 

q,k 



A k + Al k 



B k + B 



a„ +a 



a q + a'_ q 



(8) 



where the coupling constants between the photons and 
excitons take the following form 



Gi(g) = VN^gi(q) 



'N T ——±du 



VTlUJn 



G 2 (q) = VN^g 2 (q) = J N T ^j-d 2 , 



(9) 



here w q — \q\c. The wave-vector matching factor 17 in Eq. 
(8) is 



0(k + q ) = lY,e^ = l Mk f + Na \ (10) 
y 11 N ^ N sin(^a) V ' 



which is real and equal to 1 for k+q = 0, and 0(k+q) < 1, 
for k + q ^ 0. 



III. EQUATIONS OF MOTION WITH ROTATING 
WAVE APPROXIMATION 



In this section, we consider the rotating wave approx- 
imation (RWA) to deal with the interaction • A be- 
tween the radiation field and the multi-atom system. Af- 
ter introducing the two-mode exciton operators in Eq. 
(6), the interaction Hamiltonian between the excitons 
and photons with RWA is obtained as 



H RWA = hJ2G 1 (q)0(k + q) 



q,k 



A k al q + Al k a q 



+h^G 2 (q)0(k + q) 

q,k 



B k a 



(11) 



where we have neglected the higher frequency (non- 
resonant) terms: A^_ k a}_ ql A k a q , &_ k a}_ q: and B k a q . 

The Heisenberg equations for the exciton and photon 
operators are 



9 A 

d 



fiiifc + J2 G i(q)0(q - k)a q , 



i—B k = fl 2 B k + G 2 (q)0(q - k)a q , 



dt 
d 



,—a q = uj q a q + Gi (q) 0(k - q)A k 



+G 2 (q)J20(k-q)B k . 



Taking the following transform 

A k - A k = A k e l ^\ B k ^B k = B k ^\ 



(12) 
(13) 

(14) 

(15) 



to remove the fast varying factors, we obtain the formal 
equation for the exciton operator 

^A k {t) = -i^Gi(q)0(q - k)e-^-^a q {Q) 

-J2Gi(q)G 2 (q)0(q-k)0(k' - q) 

q,k' 



[ B k > (f') e -iK-ni)* e iK-n2)t' dt ' 
Jo 

Y,G 2 MO{q-k)0{k' -q) 



q,k' 
ft 



f Avit^e-^o-^V-^dt'. 
Jo 



(16) 



Here, the first term proportional to a q (0) is the so-called 
quantum noise term. The second term in the above equa- 
tion corresponds to a multi-photon processes (MPP) in- 
cluding stimulated emission and absorption effects. Its 
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contribution can be ignored since it is a higher order term 
from the standpoint of perturbation 40 . The last term in 
Eq. (16) can be solved by using the Weisskopf-Wigner 
approximation (WWA), i.e., assuming that Ak>(t') varies 
sufficiently slowly so that it can be factorized outside 
the integral. The remaining part of the time integral of 
the last term in Eq. (16) can be evaluated and get a 
Dirac 5 function with variable {u) q — f2i) and a principal 



part V 



term, which contributes a frequency shift 



-fii 

(Lamb shift). 

The equation of motion for (t) can also be obtained 
in the similar way. By using the WWA, and neglecting 
the MMP, we can solve the equations for both of the two- 
mode excitons to obtain Ak(t) and Bk(t). Substituting 
them into Eq. (14), the explicit expression for a q {t) in 
terms of a q (0), Ak(0) and Bk(0) can be obtained. Finally 
we can get the positive- frequency part of the electric field 



where we have taken the photon normalization volume 
V to be LA, with A being the area of the crystal slab, 
and assumed that the slab is located at the middle of the 
volume. When L is sufficiently large, the sum over q has 
been replaced by an integral: ^2 q Jf^, dq.... 

For an arbitrary initial state \<j>) = \<j) ex ) ® \4>l) of the 
total system, the light intensity radiated from the two- 
mode exciton system is defined as 



I(z,t) = — (4>\ : &-\z,t)E(+\z,t) 



(18) 



where the symbol ": ... :" means the normal product 
according to the exciton operators, z is the position of 
a detector relative to the crystal slab. The first-order 
degree of coherence of fluorescence is given by 

ff(1)(z;M + r) = £(*\:M-Hz,t)El+H*,t + T):\<l>) 



y/I{z,t)I(z,t + T) 



(19) 



7>\ ^fv'^'- ^ 



Similarly, one can also define the second-order degree of 
coherence as 



g^(z;t,t + r) 



in 2 



: E(~){z, t)E(~)(z, t + t)M+\z, t + T )M+) (z, t) 



I(z,t)I{z,t + r) 



(20) 



IV. QUANTUM INTERFERENCE EFFECT IN LIGHT INTENSITY 



In this section, the features of exciton fluorescence in 
the case of monolayer will be studied in detail. From the 
definition of the wave vector of excitons, when N = 1, it 
takes only one value k = and the wave- vector matching 
factor 0(q -k) = 0(k' - q) = 1. Thus with the help of 
WWA, we get the equations of motion for exciton modes 
A and B as 

|l (t) = -i^d^e-^-nO'a,^) - |loW, (21) 
^Bo(t) = -i^G 2 (g)e-^-^)*a g (0) - |s (t), (22) 

Q 

in which, 

^^^ = 3 U^ 72 ' (24) 

where, 71 = 4 3 ^ 1 is the radiative decay rate of an iso- 

lated atom from |ei) to \g), and 72 = 3 ^ 2 e3 2 is that of fa) 
to \g) , respectively. \j — for j = 1,2, denote the 



reduced wavelengthes. We find that, because of the im- 
plementation of WWA and the ignorance of multi- photon 
scattering processes (MPP) in the derivation of Eqs. (21) 
and (22), the two- mode excitons decay independently in 
an exponential rule. Even though our treatment seems to 
be bold, we can get a useful result that the decay rates of 
excitons in two dimensional crystal slab are proportional 
to the enhanced factor (A/a) 2 . In section V, we will con- 
sider a more complex case, in which a non-perturbation 
approach is used to restudy the radiative decay rates of 
the two-mode excitons without WWA, and the effects of 
MPP will be also investigated there. 

Substituting the solutions of Eqs. (21) and (22) into 
Eq. (14), we get the explicit expression of the photon 
annihilation operator a q (t) for the monolayer case: 

a q {t) =u q {t)A {0)+v q (t)B {0) 

+e-^*a 9 (0) + ^^(i)a 9 '(0), (25) 

where 

p -iu> q t _ p -r)it/2 p -iQ 01 t 

«*(*) - ~ o +? „/2 ' (26) 

UJ q — ll l+ lT]l/ 1 
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v q (t)=G 2 (q)- 



-ILUnt 



g— r)2t/2g— ifio2i 



LOq - f2 2 + %/2 



1 /"* / Z" 00 / 

(27) F A (z,t) = - dt'er^ / 2 e- inoit dw q 2cos[w q - 
n Jo Jo v ' 



«>«,«'(*) = 



Gi(g)Gi(gQe-^ f 
ay - O i +ir]i/2 



e i(ui q -Qoi)t e -r]it/2 _ ^ 

w g - O i + %/2 
same with 1 — » 2, (28) 



in which fioj = + Aj, for j = 1, 2, are the renormal- 
ized physical frequencies, and Aj are the lamb frequency 
shifts 



a, = -v r 

Jo 



dio, 



pM\gm\' 



L0 o 



(29) 



Here '"P" denotes for the Cauthy principle part, and 
p{uj q ) is the density of states of the radiation field. 

Therefore, one can obtain the positive frequency part 
of the electric field operator 



AAc 



&+\z,t) = E { +\z,t) 
+ 

where 



AAc 



F A (z,t)A (0) 



F B {z,t)B {0), 



(30) 



2-kTilo, 



1 [a,(o) £ 



q' 



(31) 



The first term in (z, t) is just the free varying photon 
field. The second term is proportional to the square of 
coupling constants. As mentioned in Eq. (17), the sum 
over the wave vector q can be replaced by an integral, 
thus two time-dependent functions F A (z, t) and Fb{z, t) 
in Eq. (30) are 



F A (z,t) = - 

n Jo 



0^2 cos 



UJ q - UJi 



; (32) 



i f°° ( z\ 

Fs(z,t) = — / dbJql COS I U!q — ) 

71" Jo V c / 



LOq - LU 2 



- (33) 



Here lo\ = f2oi — *?7i/2, 1^2 — ^02 — w?2/2. In the Ref. 17 , 
similar expressions were presented. However, the authors 
of Ref. 17 did not calculate the integrals over dui q , which 
may result in the light intensity going into the space- 
like region. In this paper, we will further calculate the 
integral over dLo q 40 ' 43 to give an explicit expression of 



xe .v-.,-.. U i A > ~-j m (34) 
Replacing the integral variable u) q by u> q — Qoii we have 



o i{u> q -n 01 )(t' -t) 



F A (z 



1 /"* 
,*) = -/ dt' 
k Jo 



e J7it'/ 2 g-* fi oit 



-fin 



x2cos 



(w g + O i) z 



,(*'-*) 



(35) 



We assume that O01 is much larger than all other quan- 
tities of the dimension of frequency 23 , so the lower limit 
of the integration can be extended to —00. Therefore, 
the above integral over dio q can be approximated by two 
S functions with variables (t' — t + ^) and (f — t — |). 
In the region z > outside the layer, performing the 
integral over t', we get 



F A (z, t) = 2e-^*-i)/ 2 e - iQ ° 1 (*-i)0 (t - 



2e - iwi( t-f )@ ( t 



(36) 



where 9 is a Heaviside function. Similarly we can also 
obtain 



F B (z,t) = 2e~ l ^ t ~^Q (t 



(37) 



From Eqs. (36) and (37), one can find that there is no ad- 
vanced propagator proportional to e _l (""^* + S) appeared 
in the positive frequency part of the electric field oper- 
ator E( + \z,t). So that the system is causal and well- 
behaved 40 . 

It should be pointed out that F A (z,t) and Fs(z,t) in 
Eq. (36) and Eq. (37) show that the electric field of 
the fluorescence emitted from the two-mode excitons is a 
temporal damped plane wave, but does not decay with z. 
A realistic electric field, however, should decrease with z, 
which happens, especially when a light propagates from 
a medium with low refractive index to that of higher one. 
Besides, a realistic electric field should not simply be a 
plane wave. However, the theoretical calculations with 
this subtle consideration will be very complicated and do 
not provide us with necessary physics in general. Our 
ideal assumption of a plane wave propagating in z may 
work well for longer wave length radiation near the crys- 
tal. 

The vacuum state of the exciton may be defined as 



4,(0) |0> = Bo(0) |0) = 0, 



(38) 



which represents that there is no exciton contained in 
the crystal slab. Then number states for the two-mode 
excitons are 



\n) A = J= [4(0) 



\m)j 



10} 



<m\ 



Bl(0) |0> . 



(39) 
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The coherent states can be formally defined as eigenvec- 
tors of exciton annihilation operators 



A (Q)\a) A =a\a) A , 



(40) 



It is noticed that this definition can work well only in the 
low excitation case. This is because the generic expansion 
of a coherent state concerns the Fock states with higher 
excitations. However, the coherent state with smaller av- 
erage exciton numbers can still approximately describe 
the quantum coherence natures of the exciton systems. 
Actually, we can understand the Fock state or the co- 
herent state in terms of the single atomic states (for the 
details please see the appendix A). In this sense, the ini- 
tial state of Frenkel exciton can be written according to 
the initial preparations of the single atomic states. A 
simplest illustration is that the vacuum state of Frenkel 
exciton is just the state formed by all atoms in the ground 
state. According to the appendix A, the Fock state of 
Frenkel exciton is a symmetrized many-atom state with 



certain atoms in the excited state. Especially, the coher- 
ent state of Frenkel exciton is an atomic ( SU(2)) coherent 
state. 

The fluorescence of the two-mode excitons will be stud- 
ied in the following of this section. If one take the initial 
state of total system as \<j>) — \4> ex ) ® |{0}), i.e., all modes 
of the light field are initially in vacuum state, EQ + \z,t) 
in the positive frequency part of the field operator of Eq. 
(30) gives zero contribution. Therefore, in the following 
calculations, we can neglect safely the "quantum noise" , 
the terms proportional to a 9 (0), and investigate only the 
influence of various initial exciton states on the fluores- 
cence of the low density excitations in the crystal slab. 

We find that the light field emitted from the two-mode 
exciton system contains two eigenmodes: f2 i and $l 02 . 
These two modes of the light field can give an interfer- 
ence if there are non-diagonal terms in the light intensity. 
To see this, we need calculate the light intensity firstly. 
Substituting Eq. (30) and its Hermitian conjugate into 
Eq. (18), the light intensity is 



I(z,t) = \F A (z,t)\ 2 (iJ(O)io(O)) + l -h^ \F B (z,t)\ 2 (bJ(O)Bo(O) 



FX(z,t)F B (z,t) (4(0)£ (0)\ + C.c. 



(41) 



where (...) denotes ensemble average, SI3 = y/Clii^, and 
V3 = \ZViV2- The first two terms in Eq. (41) represent the 
average intensities of the two eigenmodes: fioi and 
respectively. The last term proportional to gM^^ [...], 
however, gives a contribution of interference of the light 
fields. 

When the exciton system is initially in a separa- 
ble state p(0) = pa ® Pb and both pa and p B can 
be diagonal in Fock representation, there arc only two 
terms in the light intensity, which represent the contri- 
butions from the two-mode excitons, respectively. There 
is no non-diagonal term in the light intensity due to 

iJ(0)B (0)\ = 0. The light intensity is 



1 



I(z,t) = ^1^- {n /A 



4 e -m(t-z/c) 
+4 X e- V2(t ~ z/c) ] Q(t-z/c), 



(42) 



where \ = Xo(m) B / (n) A , \o = ^m/^iVi, and (n) A 
((m) B ) is the initial mean number of ^4o(So)-mode ex- 
citons. Note that xo is only determined by the intrinsic 
properties of the 2D sample. For a fixed xo, X can be 
used to describe the degree of unsymmetrically excita- 
tion within the crystal slab initially. Eq. (42) is plotted 
in Fig. 2, and our results show that I(z,t) increases 
abruptly after the propagation time t — 2ir/£li then de- 
cays exponentially and no interference pattern appears 
in this case. When (m) b = 0, i.e., the £?o-mode exci- 
tons are initially in vacuum state, as shown in Fig. 2 



(the solid line), our result will go back to the two- level 
crystal atoms case 17 due to the using of RWA and ne- 
glecting MPP in the derivation of Eq. (21) and Eq. (22), 
i.e., the mode of exciton being initially in the vacuum 
state gives no contribution to the dynamics of the total 
system. Besides, from the dot line ((to) b / (n) A = 1) 
and dashed-dot line ((to) b / {n) A = 5) of Fig. 2, we find 
that, with the increase of x, the amplitude of the light 
intensity becomes higher, which shows clearly that the 
amplitude of the light intensity is the sum of the contri- 
butions of the two eigenmodes of the radiation emitted 
from the two-mode excitons. 




G 



FIG. 2. Time evolution of light intensity I(z,t) at point 
2 = 2ttc/Q.\ for the case that the density matrix of the initial 
exciton state is diagonal in Fock representation. The solid line 
(Down): (m) B = 0; The dot line (Middle): (m) B = (n) A ; 
The dashed-dot line (Up): (m)s = 5(n)A- The units of 
I(z,t) is \(n) aTiQ.^ . t is in units of r/i/2fii = 0.01, 

^2 = 0.5f2i, 772 = iji- 



However if pa and ps are non-diagonal in Fock rep- 
resentation, e.g., the exciton system initially being pre- 
pared in a factorized coherent states \<f> ex ) = \ a ) a®\0) b> 
the light intensity becomes 



I(z, t) = -ftfii^ (n) A 0(t - z/c) Ue-^ 1 -*!^ + 4v e -" 2(t - z/c) 
8 A V 

+8Vxcos[(r! 01 - n 02 )(t - z/c) + <p] e -^+V2)(t-z/c)/2 j _ (43) 



where is the phase difference between the two coherent 
states. We find that the last term ^hfls^- {...} in Eq. 
(43) gives a nonzero contribution. The detector at z may 
register an fluorescence signal oscillating at frequency 
|Ooi — r2o2 1 - The temporal interference phenomenon may 
be observed by using a broadband detector 40 because 
with which there is no way to know the photon received 
by the detector is emitted from which modes of excitons. 
This unknowing for the "which- way" detection will result 
in the so called quantum beat phenomenon. 

Form Eq. (43), one can easily find that I(z,t) is os- 
cillating with the time evolution. The detector at z will 
receive the first peak after the propagation time t = z/c 
(In fig. 3, we take z — 2-kc/Q.i. It is not a necessary 
choice.). The influences of x on the time evolution of the 
light intensity are investigated in Fig. 3, where we have 
used <j> = 0, i.e., zero phase difference between the two 
coherent states. Our results show that if one of the two- 
mode excitons is initially in vacuum state (say (3 = 0), we 
recover the normal exponential decay of the light inten- 
sity. Besides, the amplitudes of the light intensity become 
larger with the increase of \. The result that the light 
intensity at initial stage does not go down to zero (the 
dot line and the dashed-dot line) is the consequence of 
unsymmetrically excitation, i.e., the initially generation 
of the two-mode excitons with different mean numbers. 



40 



20 



1 ! 



1 ' 



In 1 " M . 
H ' 1 1 !..! \ 
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FIG. 3. Time evolution of light intensity for the case that 
the excitons are initially in a factorized coherent state. The 
Solid line: |/3| 2 = |, 
Dashed-dot line: |/3| 
same with Figure 2. 



The Dot line: \(3\ 2 = 5|q| 2 ; The 
10|a| 2 . Other parameters are the 



All our mentioned above are separable state cases, in 
which the density matrix of the two-mode exciton sys- 
tem can be factorized initially. On the other hand, 
if the system is initially in an entangled state, e.g., 
\(j>ex) = (\0) A ® + I 1 ) a ® |0) B ), one may also ob- 
serve the quantum beat phenomenon. 

The first-order degree of coherence of the light field 
|<7^(t)| as a function of the time delay r is calculated 
by substituting Eq. (30) and its Hcrmitian conjugate 
into Eq. (19). We find that if one of the two-mode ex- 
citons is initially in vacuum state, e.g., (m) B = 0, then 
|<7^(t)| = 1 with regardless of the another mode state, 
which implies that single-mode exciton in the 2D crystal 
slab emits a complete coherent light. Mathematically, the 
complete coherent light is obtained when |.g^( T )l = 1- 
This happens when 41 

(E^-\z,t)E^(z,t + T)) oc£*(z,t)£(z,t + T), (44) 

where £(z,t) = (E^ + \z,t)), i.e., two-time correlation 
function of the total electric field may be factorized as 
the ensemble averages of the electric fields. As a special 
case, when there is only one mode exciton being excited 
initially, one can easily find that the factorized condition 
is satisfied. The physical meaning of g^(r) can be un- 
derstood by considering the visibility of the interference 
fringes, which is proportional to |<7^(t)|. A maximum 
visibility of the fringes is obtained when ^^(t)! = 1. 

If both of the two-mode excitons arc populated initially 
with p ex = pa ® Pb and both pa and ps being diago- 
nal in Fock representation, the factorized condition of 
the two-time correlation function is broken, which leads 
to < |<7^(t)| < 1. The magnitude of the first-order 
degree of coherence for this case is 



g^(z,t;z,t + T) 



yjl + 2xe-( I '2- I 'i) T / 2 cos(n 01 - rtmjr + x 2 e-^2-vi)r 
yj\ + x + xe~^ m ~ m ^ T + x 2 e _ ( , ' 2_r ' 1 ) T 



(45) 
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We find that | (r) | oscillates regularly with r, as shown 
in Fig. 4(a). The magnitude of the first-order degree of 
coherence also depends on the larger the degree of 
unsymmetrically excitation is, the smaller the oscillating 
amplitude of |<7^(t)| is. When \ - * °°j lff^( T )l tends 
to 1 . For the case of the excitons being initially in a fac- 
torized coherent state, as shown in Fig. 4(b), ^^(r)! is 
equal to 1 due to the satisfaction of the factorized condi- 
tion. 




20 40 60 80 100 120 

Time Delay 



0.6 



0) n 4 



0.2 



20 40 60 80 100 120 

Time Delay 

FIG. 4. (r) I as a function of time delay r. z = 2-kc/Q,\ , 
t = 27r/fii, ??i/2fii = 0.01, Q 2 = 0.5fii, rj 2 = r?i. (a) The case 
that the density matrix of initial exciton state is diagonal in 
Fock representation. The solid line (up) is (m)s = 50{ii)a; 
The dot line (middle) is (m)s = 10(ti)a; and the dashed-dot 
line (down) is (m)s = {ti)a- (b) For a factorized coherent 
state case. 

The second degree of coherence |</ 2 H r )l as a func- 
tion of r is also calculated. We find that, when one 
of the two-mode excitons is initially in a vacuum state 
(say (m) B = 0), then: (1) |.g (2) (r)| = 2, for the case 
of the other mode exciton being in a chaotic state; (2) 
|c/( 2 )(t)| = 1 — 1/ n, for the case of the A -mode excitons 
being initially in a Fock state; (3) |<7^(t)| = 1, for the 
coherent state case. If both of the two-mode excitons are 
populated initially, we find that |g^( T )l oscillates reg- 
ularly with r, and the oscillating amplitude of |<7^(r)| 
becomes smaller with the increase of \- 



V. NON-PERTURBATION APPROACH 
WITHOUT WWA 

In general, the rotating-wave approximation and 
Weisskopf-Wigner approximation are used frequently in 



quantum optics to study the fluorescence emitted from a 
single atom system. However, for the case of exciton radi- 
ance from a multi-atom system, we do not know exactly if 
they are valid because of the super-radiance feature of ex- 
citons. In this section, we will restudy the fluorescence of 
Frenkel excitons in the V-typed crystal slab by using non- 
perturbation approach 3 ^ 6 without the rotating-wave ap- 
proximation in the interaction Hamiltonian (e/mc)P • A. 
The self-interaction term (e 2 /2mc 2 ) A 2 is also included in 
our model to avoid un-physical roots in the characteristic 
equations for the exciton dispersion relation 17 . Both the 
stimulated emission and reabsorption effects are taken 
into account in the new theoretical treatment. We focus 
our attention on the radiative decay rates of the excitons 
and the light intensities for various exciton initial states 
to compare with our previous results in Sec. IV. 

The interaction Hamiltonian H\ between the Frenkel 
excitons and photons without RWA has been given in 
Eq. (8) in section II. Our further discussion will also 

-A into account. In 



take the self- interaction term 
second quantization, this self-interaction Hamiltonian is 
written as 



H 2 = hJ2f(q,q') 

g,g';l 



it 



a q > + al q , 



J(q+q')la 



in which 



,, M N L ire 2 ^ 
/(«>?) = ~ 7j=7f 2^ <V ' e q' v , 



(46) 



(47) 



where e qA is the unit polarization vector of the (q, X)- 
mode photon with A = 1,2. The sum J] A A , e q A • e q /,v 
had been given for the two-level atom case 44 ' 17 with cer- 
tain assumptions. Following them we can give a similar 
result for the V-type three-level case (the detailed calcu- 
lations are presented in Appendix B): 



/W) = 



2ttN l 



chVyf\qq 
1 



= r^id 2 



n 2 d 2 2 ] 



N 



l-G 1 (q)G 1 (q') + l-G 2 (q)G 2 (q') 



(48) 



By making use of 

J2 e t{q+q '^ a = NJ2 0(q' - k)0(k + q), (49) 



we can get the self- interaction Hamiltonian of photon H 2 . 
The Heisenberg equations for exciton and photon oper- 
ators can be obtained from the total interaction Hamil- 
tonian H int = Hi + H 2 (Eq. (CI) in Appendix C). Per- 
forming half-side Fourier transformation (HSFT) 



u>) = f 
Jo 



A q {t)e luJt dt 



(50) 
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on both sides of the equations, four algebraic equations 
about Ak{uS), Bk(u>), a q {uS), and a}_ q {uo) are obtained as 
listed in Appendix C. Combining these equations we get 



Eq. (51). Then we get the electric field by the inverse 
HSFT 



a q {u) - al q (uj) 
fti(w 2 -<7 2 c 2 ) 



Gi(?)^0(fc-(|) 



A k (u)-Al k (w)] -i[A k (0)-Al k (0) 

G 2 (q)J20(k-q) 



tl 2 {0J 2 - q 2 C 2 ) 
: u [B fc (w) - Bl k (u)] - i [B fc (0) - Sl fc (0) 

l TT a ^ - "rETT -^ )' 

ijj - \q\c uj + \q\c 



E(z,t) = 

where 

E(z,u) = 

= i 



1 

2^ 



oo+ie 



E(z,u)e- iuJt dio, 



(52) 



oo+ie 



2ir\q\ch 



V 



q\chL 
2irA 



a q {w) - al,(w)] 

o,(w) - al,(w)l e iqz dq, (53) 



(51) 



for the electric field of the fluorescence. In order to 
calculate it, one need to know A k {uj) — A^_ k (w) and 

Bk(co)—B^_ k (w), which are determined by a q (u))+at_ q (ui). 

In appendix C, two coupled equations for A k (ui) — A^_ k (u>) 

and Bk(uj) — B_ k (uj) are presented in detail. We solve 
these two equations consistently and replace them into 



for sufficiently large L. 

From the explicit form of the electric field operator, the 
radiative decay rates of Frenkel excitons can be solved. 
For the sake of simplicity, we will restrict ourselves to the 
monolayer case N = 1, in which only zero wave vector 
of excitons are involved in our model. Solving the two 
coupled equations for Aq(uj) — Aq(u>) and Bq(lo) — B\{uj) 
(see Eqs. (C13) and (C14) in Appendix C), we get two 
independent equations 



uj 2 - n 2 - 



2iJ 2 



lo a - Qi - 



2u^ 
2 



^ ( o 2) H 



fi 2 °° 



(w + ni)A)(o)-(w-fii)4(o) 



+2iojF$(lu) Uu + n 2 )B (0) - (w - n 2 )Bj(0) 

2 "-^H (^-n 2 ) k(o)-4(o) 



and 



(54) 



O 2 2uj2 F {1) d,\ 



1U, S ■ Bit,) 



n 2 



Q 1 



{u + n 2 )B (o) - (w-n 2 )Bj(o) 



+2i^>) [(w + Oi)Ao(O) - (w - fii)4(0) 

-i^i&VX" 2 - n?) k(o) - bI(o) 

\l 2 L 



(55) 



where 

^ ) M = -^^ = ^/^ = ^. (56) 

for j = 1, 2, 3. Fqo^(w) for j = 1, 2 give the exciton 
wave function's overlap for the A - and i?o-mode exci- 



ton, respectively, however Fq (u>) is that between the 
two modes. Note that MPP has been taken into account 
in Eqs. (54) and (55). Then, we get the same character- 
istic equations for decay rates and frequency shifts of the 
two exciton modes A and Bo 
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CM - o, 



(57) 



where 



((lj) = (lj 2 — Q 2 + iriiLjj)(uj 2 - Sl% + iri 2 Lu) + rji^LJ 2 . 

(58) 

The roots of Eq. (57) can be solved exactly, 

u)\ = Qoi - ill/2, lj 2 = — fioi — 

w 3 = «02 - ir 2 /2, w 4 - -^02 - %T 2 /2, (59) 

which determine the poles of A$(lj) — Aq(lj) and Bq(lj) — 
Bq(u>). The imaginary parts of the roots are the decay 
rates of the excitons, and the real parts are the renor- 
malized physical frequencies. For the case Oi = Q 2 = 
i.e., the degenerate case, we get 



wi = O, 0J3 = Cto - i (r/i + 772) /2, 
lo 2 = -ft, Ui = -Ct -i{vi + 772) /2, 



where 

n = \ln- 



1 2 



(m + mY 

8ft 2 



(60) 



(61) 



For the non-degenerate case fl\ ^ CI2, we assume that 
the physical roots of the characteristic equation are not 
far away from iClj due to the condition r/j -C Clj. There- 
fore, one can expand lj up to third-order of r]j/Slj, and 
get 



f2oi = ^1 



1 - 



02 





1 mv2 






2 ^2 




1 V1V2 






2 ^ 1 




-r?iO| 

„.,2 % 





r 2 = 



'/2 



r? 2 Q 2 



(ft| 



^2 ^ 



(62) 



where 77^ with j = 1,2 are given in Eqs. (23) and (24). 

We compare our results with the two-level lattice 
atoms case 17 , in which the radiative decay rate of ex- 
citons is 77 = 3 7j where 7 is the decay rate of 
an isolated lattice atom, A = c/fi denotes the reduced 



1 - 



2 

8S2 2 



wavelength, and physical frequency is Sl = SI 

However, we see from Eq. (62), for the three- level lat- 
tice atom case, the decay rates Tj are different from the 
two-level atom case rjj by a third-order correction, which 
comes from MPP (including stimulation emission, ab- 
sorption and other high-order processes). Here fioi — ^1 
and O02 — ^2 denote the frequency shifts for the ^4o-mode 
and £?o-mode exciton, respectively. For the two-level case 



the frequency shift is gjp-. Whereas, for three-level case 
considered here, the phase shifts for the two-mode exci- 

1 vim i 1 

2 n 2 -n 2 ' ' J 



ton are ^ 



1,2. The last term in the 

frequency shifts is a second-order correction, which can 
be adjusted by tuning the energy spacing of the upper 
two levels of the lattice atoms. 

In the remainder of this paper, we will present the ex- 
plicit expressions of the electric field operator. The time 
evolution and spatial distribution of the light intensity 
are calculated to compare with the results of KWA in 
Sec. IV. 

From Eqs. (54) and (55), the explicit expressions for 
A (uj) - A\(lj) (Eq. (C15)) and B (u) - Sj(w) (Eq. 
(C16)) are obtained. Replacing them into Eq. (51), we 
get the explicit light field operator in terms of the initial 
exciton field operators 



+ 



a q (uj) - al q (oj) 
2iuG 1 (q)(uj 2 - Ql) 

(lj 2 - q 2 c 2 )((uj) 
2iujG2(q)(uJ 2 - n 2 ) 

(lj 2 - q 2 c 2 )((co) 



(w + fii)4>(0) + (w-fti)4(0) 

(w + ft 2 ).Bo(0) + (u;-ft 2 ).Bj(0)] , 

(63) 



where we have neglected the quantum noise terms pro- 
portional to a q (0) and a_ q (Q) in the above calculation. 

Substituting a q (uj) — a_ q (oj) into Eq. (53) and choosing a 
proper contour integration in the upper complex q plane, 
we get E(z, lj) in the positive z region outside the crystal 
slab as the following, 

E(z,uj) 



= «' cA CH 2 ^ + ^o(0) 
+(w-ni)4(0)] e^ z 

+ 'V cA ^M ^ + n^oW 

+(w-n 2 )Bj(0)] e^ z , 

(64) 

where A is the area of the layer. As mentioned above, we 
choose the normalization volume for the photon as AL. 

The electric field E(z, t) in the positive z region outside 
the crystal slab can be calculated by Eq. (52). For the 
case z — ct > (< 0), the contour integral over du can be 
performed by choosing the integration path in the upper 
(lower) complex lj plane. The result is that E(z,t) = 0, 
for z — ct > 0; but for z — ct < 0, 



E(z,t) =£( + \z,t)+K.c. 



(65) 



The explicit expression for £^ +S) (z,t) in Appendix C 
shows that, in positive z region, the electric field gen- 
erated by the two-mode exciton is damped exponentially 



10 



with two eigen-decay rates Ti and r 2 . The corresponding 
eigen-modes are linear combinations of Aq- and i?o-mode. 
The electric field in negative z region can also be derived 
similarly and shows that it is a damped wave propagat- 
ing in the negative z direction. Since all the roots of the 
characteristic equation are in lower half plane of complex 
cj plane, then E(z, t) = 0, for z — ct > 0, so our result is 
reasonable and obeys the causal rule. 

Similarly we can also calculate the positive frequency 
part of the electric field operator, which is determined by 
a q (uj). The explicit expression of a q (u) can be obtained 
by substituting Eqs. (C15) and (C16) into Eq. (C4) 
(sec details in Appendix C.) and omitting quantum noise 
terms. We get 



lchn im (lo 2 -Q 2 ) 



1 V At: A C(w 

x 4>(o) + (w-ni)4(o) 

chVL 2 rj 2 (u 2 - ft? 



[(w + fii) 



+ 



[(w + ft 2 ) 



AnA CM 
x B (0) + (u-n 2 )Bl (0)]} 



(67) 



A straightforward calculations give the following re- 
sults: 



E {+) (z,t) = 0, for t < 0, 



(68) 



(u> - \q\c)a g (u) 



iGi(g)- 



CM 



+ iG 2 (q) 



ft 2 



CM 



(w + ni)i4o(0) + (w-fii)4(0) 
(w + ft 2 ) S (0) + (w-ft 2 ) Sj(0) 



E {+) ( Zl t) 



AAc 

?(-) 



F { A + \z,t)A o (0) 



+F { -\z,t)Al(0) 



Therefore, we get 



dq 



oo-\-ie 

dco 

-oo+«e 

oo+ze 



q\chL 
~2t7a~°" 1 

„-iwt 



a q (u))e 



(66) 



i(qz-ujt) 



nUCl 2 r] 2 



AAc 

+F£\z,t)Bl(Q) 



F B +) (z,t)B o (0) 

, for t > 0, 



(69) 



duj- 



oo+ie 



U>„ 



■ UJ 



in which the time-dependent coefficients F^ z \z,t) and 
F B ±] (z,t) are 



/•GO 

i / 


dio q 2 




-lUJqt 




t Jo 


((Jg - Wl) (Wg - W 2 ) (Wg - LU 3 ) (uj q 


- w 4 ) 




H 


-ft 2 .) (wi ± fii) e"^ 1 * 






H 


- uji) (u>i - lo 2 ) (wi - w 3 ) (wi - w 4 ) 










-ft 2 ,) (w 2 ±0i)e-^ 2t 






K 


- W 2 ) (<*>2 — Wl) (U! 2 — UJ 3 ) (uj 2 — OJ4) 










-ft 2 ,) (w 3 ±f7i) e -^ 3t 






K 


- u 3 ) (lu 3 - U!l) (oj 3 - u> 2 ) (lo 3 - UJ4) 








H 


-ft 2 .) (w 4 ± ft x ) e"^ 4 * 


COS 


h£ 


K 


- OJ4) M — Wl) (W4 — W2) (W4 — W 3 ) 


i r°° 


div q 2 


(^ 2 - ft 2 ) k ± ft 2 ) e - 


-iUJqt 




* Jo 


(tJg - Wl) (Wg - LU 2 ) (U>g ~ LU 3 ) (uj q 








-ft 2 ) (wi ± fi 2 ) e - ^ 1 * 






K 


- ui\) (u>i - lo 2 ) (ui - uz) (ui - LO4) 










-ft 2 ) (w 2 ±ft 2 )e^ 2t 






K 


— oj 2 ) (uj 2 — UJi) (uj 2 — U>z) (oJ 2 — W4) 








(-1 


- ft 2 ) (w 3 ± ft 2 ) e"^ 3t 







(70) 



W2) (w 3 - w 4 ) 

\ —i.ti)A-f. 



(u 2 - ft 2 ) (o> 4 ± ft 2 ) e-^ f 

(Uj q - UI4) (UJi - UJl) (u)4 - UJ 2 ) (0J4, - UJz) 



COS (uq^J , 



(71) 
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where uij, for j = 1, 2, 3, 4, are four roots of the charac- 
teristic equation of Eq. (57). We again need to carry out 
the integration over dui q in the above two equations. It 

is easy to prove that F A + ^ (z, t) can be simplified as 



F A + \z,t) 



dto q 2 cos(w 9 



2> 



ci- 



-iu) Q t 



+ c 2 - 



+ c 4 - 



W, - ^2 



c 

+ c 3 - 



LO q - LUl 

lUJqt 



LU„ — LU4 



(72) 



where the four coefficients are 



Cl 



C2 



C3 



(o>i - w 2 ) (wi - W3) (<*>i - w 4 ) ' 

(ujI - nfj (uj 2 + cii) 

(u>2 - UJl) {UJ2 - OJ3) (u>2 - W4) ' 



til) (lo 3 + til) 



(u> 3 - U)i) (u> 3 - C0 2 ) (W3 - W 4 ) ' 

H - "1) K + no 

C4 = 7 T7 T7 r- (73) 

(CJ4 - U)l) {U>4 - LV2) (W4 - W3) 

The above integral over c?o> 9 can be done as Sec. IV. The 
explicit expressions of F t " A \z 1 t) and F B ± \z,t) can also 
be obtained with the same calculations. Our results also 
confirm the causal rule. 

For an arbitrary exciton initial state \4> e x)i the light 
intensity is defined by Eq. (18). When the excitons are 
initially in a state with density matrix p(0) = pa® Pb, 
and both pa and ps diagonal in Fock representation, the 
light intensity is 



I(z,t) 



-hti 2 — {m) B 



F A + \z,t) 



+ 



?(-) 



(z,t) 



p(+) 



(z,t) 



?(-) 



(z,t) 



(74) 



which is shown in Fig. 5. When (7x1)3 = 0, the solid 
line, our result will go back to the two-level lattice atoms 
case 17 . From the dot line ({m) B / (n) A = 1) and dashed- 
dot line ({m) B / (ri) A = 5) of Fig. 5, we find that with 
the increase of %, the amplitude of the light intensity be- 
comes higher, which is the same with that of obtained 
by using RWA (Fig. 2). Contrary to the results of Sec. 
IV, however, when we consider it without RWA, the light 
intensity does not decay exponentially but in an irregu- 
lar way due to the existence of counter-rotating terms 
in E( + \z,t). Besides, it's also deserved to mentioned 
that the contributions of the non-rotating terms do not 
appear as quivers presented in Ref. 17 . 



FIG. 5. Time evolution of light intensity I(z,t) at point 
z — 2nc/Qi for the case that the density matrix of the ini- 
tial exciton state is diagonal in Fock representation. The solid 
line (Down): (m) B = 0; The dot line (Middle): (m) B = {n) A ; 
The dashed-dot line (Up): (m)s = 5(ri) A - Other parameters 
are the same with Fig. 2. 

For the case that the two-mode excitons are initially in 
a factorized coherent states, the light intensity becomes 




I(z,t) 



8 A 



H 2 F A + \z,t) +H 2 F { -\z,t) 



HarF A + >(z,t)FP(z,t) + a'F A ->*(z,t)F A ^>(z,t) 



A-) 



2 !?(-)*, 



F B + \z,t) 



7 B~\*,t) 



+ (n 2 F B +> (z, t)F^ (z, t) + (3 2 F B ->(z, t)F B +) (z, t) 
+\hti 3 ^ [a*f3F A +) *(z,t)F B + \z,t) + a*f3*F A +) *(z,t)F B -\z,t) 
+af3F A -> (z, t)F B +) (z, t) + af3*F A ~>(z, t)F^ (z, t) + C.c] , 



(75) 



where the last term ^ ?j,r2 3 [. . .] gives a temporal inter- 
ference term. The effects of phase difference between the 
two coherent states on the light intensity are studied in 



Fig. 6. We find that the first peak (at t = z/c) of the light 
intensity becomes lower in amplitude with the increase 
of the phase difference from <j> = to <j> = tt (monotonic 
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regime). In fact, the whole curve will be shifted left with 
the increase of the phase difference within the monotonic 
regime, which leads to the magnitude of the peak (at 
t = z/c) becomes lower (comparing Fig. 6(b) and Fig. 
6(c) with Fig. 6(a)). When <f> = n, I(z,t = z/c) tends 
to zero, i.e., the first peak is disappeared (see Fig. 6(c)). 
Our results also show that both the phase difference and 
the degree of unsymmetrically excitation do not effect 
the oscillation frequency of the light intensity which is 
determined by the exciton splitting. The oscillation be- 
havior in the light intensity may take place as long as 
|Sloi — r2o2 1 > Tj, i.e., the exciton splitting is greater 
than the natural lincwidth of exiton, so that one can ob- 
serve the beating phenomenon within the lifetime of the 
exciton. 




FIG. 6. Time evolution of light intensity for the case 
that the excitons are initially in a factorized coherent state: 
(3 — ae'*, with <j> being the phase difference between the two 
coherent states, (a) = 0; (b) <f> — it/2; (c) <f> = 7r; Other 
parameters are the same with Fig. 2. 

The spatial distribution of the light intensity is plot- 



ted in Fig. 7 for the case that the two-mode excitons 
are initially in a factorized Fock state. Our results show 
that I(z, to) increases exponentially (see Eq. (42)) within 
the regions < \z\ < ct for the case of rotating- wave 
approximation, and vanishes immediately as \z\ goes be- 
yond cto- The solid lines in Fig. 7 are obtained without 
RWA and shows small-amplitude oscillations due to the 
contribution of counter-rotating terms. Compared with 
the results of Ref. 17 , our results show that not only the 
electric field E(z, t) but also the light intensity I(z, t) do 
meet the requirement of causal rule. 
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FIG. 7. Space distribution of intensity I(z,to) with 
(a)t = 2%/fli, (b)i = 37r/fii, and (c)to = 57r/fli. 
(m)s = (n)A, z — ct is in units of c/fh. Other parameters 
are the same with Fig. 2. The dot line is obtained by using 
RWA, and the solid line is obtained without RWA. 



VI. CONCLUSION 

In summary, we have studied the collective radiations 
of a collection of many V-type three-level atoms in a crys- 
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tal slab. By introducing two-mode exciton operators in 
the large N limits of the collective quasi-spin operators, 
these coherent radiations can be depicted as the fluo- 
rescence of low density Frenkel excitons. The exciton 
fluorescence exhibits the stronger coherence natures that 
the statistical characters of spectrum are identical from 
the initial to final stages. This is indeed different from 
the ensemble situation with free-moving atoms that the 
atoms need a finitely-long time to produce a cooperative 
radiation, the enhanced fluorescence. 

As a main result of this paper, the occurrence of the 
quantum beat aroused from which quantum states of 
Frenkel exciton is investigated in this paper. Our results 
show that not all the quantum states of excitons may 
lead to the oscillating behavior in the light field radiated 
from the two-mode excitons. This quantum interference 
phenomenon may be observed when the two-mode ex- 
citons arc initially in a factorized coherent state or an 
entangled state. We expect that our theoretical study of 
quantum beat would be helpful in practical experiment 
to measure quantum states of excitons. Further study of 
the anharmonic exciton-exciton interaction in the model 
is needed. It is also pointed out that, the algebraic con- 
sideration for the definition of the multi-mode excitons 
can be generalized to study other exciton system, e.g., 
the quasi-spin wave collective excitations of A-type atom 
collection in lattice of crystal that can be used as a new 
type quantum memory. 
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APPENDIX A: SU(3) ALGEBRA STRUCTURE 
OF EXCITON OPERATORS FOR MANY-ATOM 
SYSTEM 



We consider an ensemble of N two-level atoms with 
their ground states \g)i and the excited ones |e) J . Since 
we can define the quasi-spin with the Pauli operators 



«r = ifl>"H 



t(e\-\gf(g\, (Al) 



the total quasi-angular momentum operators 



N N 

J_=£aL,J + = (J-V\j, = -5>;, 



(A2) 



define a representation with the highest weight J = N/2. 
The (2 J + l)-dimensional irreducible spinor representa- 
tion of SU(2) in the symmetric subspace is embedded in 
the total Hilbert space of dimension 2 N . 

For a physics system with dynamical SU(2) symmetry, 
its Hamiltonian H = H( J_, J + , J z ) is a functional of J± 
and J z . Since the Casimir operator J 2 commutes with 
J± and J z , the eigenvalue J (J + 1) will keep conserva- 
tion in the time evolution. Here, J can take one of the 
integers and half-integers ^, ^ — 1, ...,0. In a physical 
process, which one of these J takes depends on the ini- 
tial state of the atomic ensemble. The symmetric state 
| J = f ,M = -f ) = nili W represents the "conden- 
sate" with all atoms filling in the ground state. It is very 
similar to the situation of an electronic system that the 
filling in ground state forms the Fermi surface. In this 
sense, we can introduce the atomic collective excitation 
operators 



j- 



(A3) 



which is very similar to the exciton operators for an 
electron-hole pair. Considering that 



1 



sJ(N 



+ l) 2 - 4J+J_ 



(A4) 



as one solution for the basic angular momentum relation 

* (A5) 



Ji + J', + J; ' = T I T + 1 I , 

J+j- = Jx + Jy + Jz = "77" (77- 



J 2 Z +J Z , (A6) 



From a point of view based on the representation of Lie 
algebra, this section describes the mathematical origin of 
definition of exciton operator. Physically, this descrip- 
tion will clarify why the conception of the exciton based 
on the collective operators of many atoms can be valid 
only in the case of low excitation. We first discuss SU(2) 
algebra structure of excitonic operators for the two-level 
many-atom system since the SU(3) case for the 2D crys- 
tal slab containing V-type three-level atoms shares the 
same basic idea as that of SU(2). The detailed discus- 
sions for the SU(3) case will follow that of SU(2) in this 
appendix. 



the commutation relation for atomic collective excitation 
operators 



B, B 



(A7) 



1 for 



will be reduced to the bosonic relation 

those angular momentum state | J, M) with much smaller 
M in comparison to N = 2 J. That is to say, in the limit 
with very large N and low excitation, the collective ex- 
citation behave as a boson and thus we call it atomic 
exciton. The detailed proof for this was given in rcf. 24 
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by considering the physical realization of q-deformed bo- 
son algebra for a very large, but finite N. The main 
point to prove that is to consider that 



Jz 

N 




(A8) 



approaches —1/2 for the infinite N and the low excita- 
tion. Since the operators can make sense by acting on 
the symmetric space, only for those low excitation states 
\J,M) with very small M can let the value of J+J-/N 2 
approach zero, so that the bosonic commutation relation 
is obtained. We can also prove that 



JV 



&B = J2\ef(e\+0(- 

i=i v 



(A9) 



which means that the free part of the many-atom Hamil- 
tonian can rationally be described as a free boson in the 
large N and the low excitation. 

Because the "condensate" in a ground state plays the 
crucial role in defining the exciton operators, the intro- 
duction of exciton operators has to depend on the config- 
uration of the atoms and thus on the form of interaction. 
This is just the line, along which we will define the ex- 
citon operators for the V-type atomic system with the 
interaction Hamiltonian 39 



N 



N 



Hj = 9l tf \9f(ei\ + 92^ l^'M + H c - ( A1 °) 



i=i i=i 
It is easy to prove that 

N N 



E 1 = Y,\^Y l (9\,F l =Y,\9T^ 

t=\ i=l 
N N 

£2 = £isr< e2 i,F 2 = £| e2 n ff |, 

i=l i=l 
1 N 

#i = oE(i e in e iMsr<9i), 
1 i=i 

1 N 

H2 = ^(\9) ii (9\-\e2) ii {e 2 \), 



(All) 



generate a SU(3) algebra with the Cartan subalgebra 
spanned by Hi and H 2 . The basic commutation rela- 
tions are 



[H 2 ,E 2 ] 

[H 2 ,Ei] 

[Hi,E 2 ] 



E 2 , [H 2 , F 2 ] 



-Fi, 
-F 2 , 



-?±,[H 2 ,F 1 ] = 1 ^, 



(A12) 



Two sets of the operators {Ei, F\, Hi} and {E 2l F 2 
,H 2 } generate two non-commutation SU(2) subalgebras, 
respectively. Actually, the above four collective opera- 
tors define a spinor realization for the symmetric repre- 
sentation of SU(3) in Nl + 1 dimensional space. The 
atom number Nl determines the dimensions Nl + 1 of 
representations. In this sense, we can understand the 
two-mode excitation in terms of the large Nl limit of 
representations of SU(3), which just corresponds to the 
low density excitation region. Thus we can define the 
collective operators 



B = — E 2l = — F 2 . 
s/N VN 



(A13) 



For a very large N and low excitation, it is easy to prove 
that A and B commute with each other and obey the 
standard bosonic commutation relation. 

Now we can consider the construction of Frenkel exci- 
tonic states. In this way, the initial conditions for Frenkel 
exciton can be given in terms of the single atom prepara- 
tions. For example, when n = 0, 1, 2, the Fock states 
of the A-mode excitons 



i(' 4, )"'°> 



N 



N n l 2 A 



(A14) 



i=l 



take the symmetric excitation states 

|0> = \g,g,...g), 
1 N 

= -j=^2\g,g,...,e{,...,g), 

.7 = 1 



|2>. 



1 N 

—f= E \9,9,-,e{,..,e^,...,g), 



NV2 



3,k=l 



(A15) 



The second example is the coherent state of the A-mode 
Frenkel exciton 



\o)a oc exp(aA' t )|0) 

= [J [cosd\ g y +sin0e^|ei) 3 



(A16) 



H 



and a = \a\e 1 ^. The coherent nature 



where tan# 

of this many-atomic state is reflected by the fact that 
both and <p are independent of the index j of atoms. 
The quantum states of i?-mode excitons can also be con- 
structed by using the same procedure as discussed above. 
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APPENDIX B: THE SELF-INTERACTION TERM 
OF THE LIGHT FIELD 

In this section we will calculate e q \ -e^v for the three- 
level case. For an atom with a complete set of eigenvec- 
tors {|n)}, we have 

e g A • e q >\> = (g\e q \ ■ e q >\<\g) 

= ~%J2^ e i X ,x l n )( n IP' e «'A'lfl) 

n 

-^2(9\e q \ ■p|n)(n|x-e 9 ^'|fl') 

n 

2m 



he 2 



^2rt n (g\e qX ■ d|n)(n|d • e q ,y\g), (Bl) 



where we have used [x, p] = ih and p = 
(xHa — HatcJ . Ha is the free atomic Hamiltonian 

and gives the eigenvalue equation: Ha\ii) = E n \n). 

For the V-type three-level case, we take the three-level 
approximation in the above equation as ref. 17 for the two 
level atom case, obtaining 



e 9 A • e q '\> 



2m 
he* 



[Hi (e q \ ■ di) (di-e 9 'v) 



+tt 2 (e qX ■ d 2 ) (d 2 -e 9 -A')] > 



(B2) 



where il 1 = (E ei - E g )/h and fl 2 = (E e2 - E g )/% 
are atomic transition frequencies for \g) <-> |ei) and 
\g) <-> |e 2 ), respectively. d x = (ei|d|^> = (g|d|ei) and 
d 2 = (e 2 |d|(7) = (g|d|e 2 ) are the corresponding transi- 
tion dipolc moments. Choosing e q \ as the following 



e q i ■ di = di, e q i ■ d 2 
e q >2 ■ di = 0, e q < 2 ■ d 2 = 



: d2, 
0, 



(B3) 



then we get 



2m 



e q \ ■ e q >x> ~ ^[Hidi^A.i^A.A' + n 2 d 2 6x,iS\,x>]- (B4) 

Substituting J2\,\> e qA-e q 'A' into f(q, q'), we get Eq. (45) 
in section V. It is noticed that, for the case with single di- 
rection polarization of light, J2\ y e qA -e q 'A' = 1 strictly. 
The cut-off of the complete relation for the sum by only 
three levels, however, will lead to the departure from 1. 
Only with this cut-off approximation the exactly-solvable 
mode is built for the two mode excitons coupling to the 
quantized electromagnetic fields. 



APPENDIX C: DETAILED 
NON-PERTURBATION CALCULATION 

In this Appendix, we will add the necessary details and 
list the more expatiatory expressions for section V. We 
start from the total interaction Hamiltonian 



flint = h J2 Gi (q)O(k + q) [A k + Al 

q,k 



a a + a 



hY,G 2 (q)0(k + q) 



q,k 



B k + B'_ 



+ ^ E i-G^diq'W - k)Q(k + q) 



q,g' ,k 

< \a q + a)_ q a q , + a)_ q , 



+ ^ E <^G 2 (q)G 2 (q')0(q' - k)Q(k + q) 



q,q',k 



a q > + a ] _ , 



(CI) 



which includes the non-RWA terms and the self- 
interaction of the light field. 

The half-side Fourier transformation (HSFT) of the 
Heisenberg equations governed by the total Hamiltonian 
(CI) are 

(w - fii)Afc(w) = E Gi(q)0(q - k) a q {uj) + al q (u) 

+iA k (0), (C2) 
(w - n 2 )B fc (w) = J2 G2{q)0{q - k) \a q {uj) + al q {uo) 



+iB k (0), 



(C3) 



for excitons, and 



(oj - \q\c)a q (uj) 



Gi{q)^20(k-q) A k (w)-Al k (u) 



- ^Gi(q) E °( fc ~ «) [A fc (0) - Al k (0) 

k 

+ £-G 2 (q) 0(k - q) \B k {w) - Sl fc (w) 



'n 2 



G 2 {q)Y,0{k - q) [B fc (0)-Bl fc (0) 



+ ia q (0), 



(w + \q\c)a)_ Auj) 



(C4) 



UJ 



Gi(g)^0(fc-g) \A k (u) - Al k (u) 



UJ 



G 2 (q)Y,0{k - q) \B k (uj) - Bl k {uj 



+ ^G 2 {q) E 0(k q) \B k (0) Bl k (0) 



(C5) 
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for photons. Here, we have eliminated J2 q > G(q')0(q' — _|_ ^ y^ F^ 3 ?(u) ylfc/(u;) — (w 



fc) a g / (w) + a^ ? , (u) in the derivation of Eqs. (C4) and 
(C5). 

Combining the above two equations, we get 

(ui 2 - q 2 c 2 ) a q (ui) + a'_ q (w) 
= i(u + \q\c)a q (0) + i(oj - \q\c)at_ q (0) 

+2|g| C ^-G 1 (g) ]T 0(k - q) [A k (cj) Al k {u>) 

k 

-i^G 1 {q) Oik q) [A k (0) Al k (0)_ 

k 

+2\q\c^-G 2 (q) £ 0(k - q) [B k (w) - B^lo) 
2 k 

-i*G 2 («) ]T °( k «) [^(°) - B -kW] ■ (C6) 

k 

Substituting Eq. (C6) into Eq. (C2) and (C3), after a 
straightforward calculation, we obtain two coupled equa- 
tions for the Frenkel exciton operators 



(C8) 



E 



x \A k ,(w)-Al k ,(w) 



= i [(w + Oi)Afe(O) - (w - ni)Al fc (0) 



+ 2iw^Gi(g)0(g - fc) 



a,(0) + ol,(0) 



lu - \q\c lu + \q\c 



"1 fc , 



(C7) 



and 



E 



x [s fc ,(w)-Bl fc ,(o;) 
i [(w + n 2 )B k (0) - (w - fi 2 )si fc (o) 

^ [w-|g|c w+|q>|c 

2^E F ^H [^(0)-Bl fc ,(0)' 



1 it' 



where only the initial photon operators and the exci- 
ton operators are concerned. Three factors introduced 
in Eqs. (c7) and (c8) are 



(1) _ x - 2\ q \cG 2 (q)Q(q - k)Q{k' q) 

*kU\ u )-Z*, Lu 2 -q 2 c 2 
q 

(2) ^_^ 2\ q \cG 2 (q)Q( q -k)0(k' -q) 
{UJ) ~ 2f rf-fc* . 

^(a.) = j ^oc^ - g) a (C9) 

q LU q C 

which represent the overlap of exciton wave functions 
with different wave vectors. We take the photon nor- 
malization volume V to be AL where A is the area of 
the crystal slab, and place the slab at the middle of the 
volume. When L is sufficient large, the sum over q can 
be replaced to a integral: J2 q — ~* ^ S°°oo dq---- Thus 



Nana 2 

47TC 2 



2 roc 



I 



(cio, 



By carrying out the integrations as in Rcf. 17 , the ex- 
plicit expressions are 



F kk' M 



SNcw 



sin^±^7Va ^'^ N " 



sin ■ 



sin -^TVa e l ^ Na 



sin ■ 



+ 



sin K^Na 



sin - a 



sin^a S in(^±i a ) s in(^a) 



(Cll) 



2 "V 2 

When iV — > oo, the first two terms tend to zero, thus 
af 2 Qi sin 



lim 



8cw S i n (^+l a ) sin (^i a ) 



(C12) 



The above results will be used to determine a q (co) — 
a^_ q (uj) explicitly in Sec. V. 

For the single lattice layer case, ignoring the quantum 
noise terms proportional to a q (0) or a_ q (0) induced by 
the background light field, we obtain the coupled equa- 
tions of the exciton operators 



2^_ 

n 2 



= I 



(cj + n 2 )B (o) - ( w -n 2 )Bj(o) 
n 2 



-2i^-F${cu) fs (0)-S5(0) 
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-2i£-F$(u) (A)(0)- 4(0) 

+^ ( o 3) H [am -4m" 



(C13) 



and 



w — sZj_ — -^j-^oo M 



am -4m 



= 2 



(w + fii)i4o(0)-(w-fti)4(0) 



-2i^-F ( g ) H (A)(0) -4(0) 
-2»^F ( 3) H (b (0)-B5(0) 



(C14) 



These two equations lead to two decoupled equations Eq. 
(54) and Eq. (55) in Sec. V for AH - 4M and 
B (ui) - Sj(w). The solutions of Eq. (54) and Eq. (55) 
are 



AM -4M 



^ 2 -O 2 



ir) 2 uJ 



CM 



(w + Qi)A(0)- (w-Oi)4(0) 



+ 



03773 
CM 

CM 



(^ + O 2 )B (0) - (w-O 2 )4(0) 



-n 2 ) A(o)-4(o) 



(C15) 



and 



BoM-bJm 



.O) 2 - Of + Z77lW 



+ 



CM 
o 3 ~?3 r 



(co + r! 2 )Bo(0) - (w-fi 2 )Bj(0) 



CM 
CM 



(w + O0A(0) - (w-fii)4(o) 



M-O 2 ) (Bq(O)-BJ(O) 



(C16) 



where (M is defined in Eq. (58). The above two equa- 
tions determine a q (u) — aL q (cj) and give the non-zero 
electric field E(z,t) (Eq. (65) in Sec. V) 

E(z,t) =£(+)(z,t) + H.c., 
for z — ci < 0, where 

£ { +) (z,t) 

[m + 00 A(o) + M - no 4(0)" 

M - w 2 ) (wi - W3) M - W4) 

M + no A(o) + (ws - fii) 4(0) 



+ 



cA 



M 2 



2x e -i W3 (t-f) 



+ 



M - LUl) M - CJ 2 ) M - W 4 ) 

(wi + n 2 ) B (o) + (wi - n 2 ) 4(o) 



cA 

(w 2 - Of ) e -^(*-t) 



+ 



(wi - w 2 ) (wi - w 3 ) M - w 4 ) 

M + 2 ) B (0) + (wi - 2 ) 4(0) 



7r?i0 2 ?7 2 



M 



Of)< 



M - ui) M - w 2 ) M - w 4 ) ' 



(C17) 



here uij, for j = 1,2,3,4, are four roots of the charac- 
teristic equation of Eq. (57). Eq. (CI 7) shows that, in 
positive z region, the electric field generated by the two- 
mode exciton is damped exponentially with two eigende- 
cay rates Ti and T 2 . Unlike the results of Sec. IV, the 
corresponding eigenmodes, however, are two linear com- 
binations of the A) an d B modes. This is because we 
include non-RWA terms and MPP in the derivation of 
Eq. (54) and Eq. (55) in Sec. V. 



* Present Address: National Lab of Magnetism, Institute 
of Physics, Chinese Academy of Sciences, Beijing 100080, 
China. 

a Electronic address: suncp@itp.ac.cn 

b Internet www site: http:// www.itp.ac.cn/~suncp 

1 J. J. Hopfield, Phys. Rev. 112, 1555 (1958). 

2 V. M. Agranovich and O. A. Dubovsky, JETP Lett. 3, 223 
(1966). 

3 E. Hanamura, Phys. Rev. B38, 1228 (1988) 

4 T. Itoh, T. Ikehara, and Y. Iwabuchi, J. Lumin. 45, 29 
(1990). 

5 L. C. Andreani,F. Tassone, and F. Bassani, Solid State 
Commun. 77, 641 (1991). 

6 J. Knoester, Phy. Rev. Lett. 68, 654 (1992). 

7 J. Knoester, J. Lumin. 53, 1001 (1992). 

8 G. Bjork, S. Pau, J. Jacobson, and Y. Yamamoto, Phys. 
Rev. B50, 17336 (1994). 

9 G. Bjork, S. Pau, J. Jacobson, H. Cao, and Y. Yamamoto, 
Phys. Rev. B52, 17310 (1995). 

10 V. M. Agranovich, D. M. Basko, and O. A. Dubovsky, J. 
Chem. Phys. 106, 3896 (1997). 

11 S. D. Boer and D. A. Wiersma, Chem. Phys. Lett. 165, 
45(1990). 

12 H. Fidder, J. Knoster, and D. A. Wiersma, Chem. Phys. 
Lett. 171, 529 (1990). 

13 Ya. Aavisoo, Ya. Lippmaa, and T. Reinot, Opt. Spetrosc. 
62, 706 (1987). 

14 B. Devead, F. Clerot, N. Roy, K. Satzke, B. Sermage, and 
D. S. Katzer, Phys. Rev. Lett. 67, 2355 (1991). 

15 R. Bonifacio and A. L. Lugiato, Phys. Rev. A 12, 587 
(1975). 



18 



16 Y. Kaluzny, P. Goy, M. Gross, J. M. Raimond, and S. 
Haroche, Phys. Rev. Lett. 51, 1175 (1983). 

17 C. Q. Cao, H. Cao, and Y. X. Liu, Phys. Rev. B62, 16453, 
(2000). 

18 Y. Yamamoto, R. E. Slusher, Phys. Today 46, 66 (1993). 

19 K. Nishioka, K. Tanaka, I. Nakamura, Y. Lee, and M. 
Yaminishi, Appl. Phys. Lett. 63, 2944 (1993). 

20 Y. Yamamoto, S. Machida, and G. Bjork, Phys. Rev. A 44, 
657 (1991). 

21 E. Burstein and C. Weisbuch, Confined Electrons and Pho- 
tons: New Physics and Application (Plenum, New York, 
1995). 

22 C. Weibuch, M. Nishioka, A. Ishikawa, and Y. Arakawa, 
Phys. Rev. Lett. 69, 3314 (1992). 

23 Y. X. Liu, C.P.Sun, and S.X.Yu, Phys.Rev. A 63,033816 ( 
2001). 

24 Y. X. Liu, C. P. Sun, S. X. Yu, and D. L. Zhou, Phys. Rev. 
A63, 023802, (2001). 

25 Y. X. Liu, N. Imoto, K. Ozdemir, G. R. Jin, and C. P. Sun, 
Phys. Rev. A 65, 023805 (2002) 

26 M. D. Lukin, S. F. Yelin, and M. Fleischhauer, Phys. Rev. 
Lett. 84, 4232 (2000). 

27 M. Fleischhauer and M. D. Lukin, Phys. Rev. Lett. 84, 
5094 (2000). 

28 M. Fleischhauer and C. Mewes, (quant-ph/0110056). 

29 M. Fleischhauer and M.D. Lukin, (quant-ph/0106066). 

30 M. D. Lukin, M. Fleischhauer, R. Cote, L. M. Duan, D. 
Jaksch, J. I. Cirac, and P. Zoller, Phys. Rev. Lett. 87, 
037901 (2001). 

31 C. P. Sun, H. Zhan, and X. F. Liu, Phys. Rev. A58, 1810 
(1998). 



32 C. P. Sun, S. Yi and L. You, Phys. Rev. A63, in press 
(2003). 

33 Y. Li, S. Yi L. You and C. P. Sun, Science in China (series 
A), in press (2003). 

34 M. G. Raizen, R. J. Thompson, R. J. Brecha, H. J. Kimble, 
and H. J. Carmichael, Phys. Rev. Lett. 63, 240 (1989). 

35 Yifu Zhu, Daniel J. Gauthier, S. E. Morin, Qilin Wu, H. 
J. Carmichael, and T. W. Mossberg, Phys. Rev. Lett. 64, 
2499 (1990). 

36 T. Pellizzari, S. A. Gardiner, J. I. Cirac, and P. Zoller, 
Phys. Rev. Lett. 75, 3788 (1995). 

37 B. Julsgaard, A. Kozhekin, and E. S. Polzik, Nature 413, 
400 (2001). 

38 J. Knoester, Frank C. Spano, Phys. Rev. Lett. 74, 2780 
(1995). 

39 C. P. Sun, S. X. Yu, Y. B. Gao, (quant-ph/9809079); C. 
P. Sun, D. L. Zhou, S. X. Yu, Y. X. Liu, and Y. B. Gao, 
(unpublished) 

40 P. Meystre, and M. Sargent III, Elements of Quantum Op- 
tics (Springer-Verlag, Berlin Heidelberg New York, 1999). 

41 M. O. Scully, and M. S. Zubairy, Quantum optics (Cam- 
bridge University Press, Cambridge, 1997). 

42 E. Hanamura, J. Phys. Soc. Jpn. 37, 1545 (1974); ibid, 37, 
1553, (1974). 

43 W. W. Chow, M. O. Scully, and J. O. Stoner, Jr., Phys. 
Rev. A 11, 1380 (1975). 

44 P. W. Milonni, R. J. Cook, and J. R. Ackerholt, Phys. Rev. 
A40, 3764 (1989). 

45 C. P. Sun, and H. C. Fu, J. Phys. A. Gen. Math. 23, L893 
(1989). 



19 



